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ABSTRACT: For a system consisting of a large collection of particles, a set of variables 
that will generally become effectively classical are the local densities (number, momentum, 
^ . energy). That is, in the context of the decoherent histories approach to quantum theory, 
\ it is expected that histories of these variables will be approximately decoherent, and that 
^ their probabilites will be strongly peaked about hydrodynamic equations. This possibility 

is explored for the case of the diffusion of the number density of a dilute concentration 
of foreign particles in a fluid. This system has the appealing feature that the microscopic 
', dynamics of each individual foriegn particle is readily obtained and the approach to local 
■ equilibrium may be seen explicitly. It is shown that, for certain physically reasonable 
initial states, the probabilities for histories of number density are strongly peaked about 
evolution according to the diffusion equation. Decoherence of these histories is also shown 
for a class of initial states which includes non-trivial superpositions of number density. 
Histories of phase space densities are also discussed. The case of histories of number, 
momentum and energy density for more general systems, such as a dilute gas, is also 
discussed in outline. When the initial state is a local equilibrium state, it is shown that the 
histories are trivally decoherent, and that the probabilities for histories are peaked about 
hydrodynamic equations. An argument for decoherence of more general initial states is 
given. 



1. INTRODUCTION 



One of the aims of the decoherent histories approach to quantum theory is to supply an 
account of the emergence of an approximately classical world from an underlying quantum 
one [1,2,3,4]. To date, most of the effort in this direction has concentrated on systems 
in which there is a natural separation into distinguished system and environment. In the 
quantum Brownian motion model, for example the distinguished system is a large particle, 
and the environment a bath of harmonic oscillators. Histories of imprecise positions of 
the distinguished system are then typically decoherent as a result of their interaction with 
the environment, and the probabilities for histories of position are typically peaked about 
classical equations of motion, with dissipation [3,5]. These models have the advantage of 
technical simplicity - a variety of known techniques, such as path integral and the influence 
functional [6,7], may be brought to bear, and emergent classicality is readily shown. These 
systems have also been analysed from other perspectives, including the density matrix 
approach and quantum state diffusion approach, with qualitatively similar results [8,9,10]. 

It is however, of considerable interest to investigate models of a more general class that 
might describe a gas or fluid, where there are certainly no distinguished point particles, 
and possibly even no obvious candidates for an environment to produce decoherence. In 
this context it has been suggested by Gell-Mann and Hartle that for a large and possibly 
complex system, the variables that will become classical "habitually" are the local densities 
(energy, momentum and number, for example), integrated over small volumes [3]. They are 
expected to be approximately decoherent because they are approximately conserved (since 
exactly conserved quantities are exactly decoherent [11]). Furthermore, there is some hope 
that the probabilities for such histories might be peaked about interesting deterministic 
evolution equations. This is because it is known that the local densities, averaged in a 
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local equilibrium state, obey hydrodynamic equations. 

It is therefore of interest to study the decoherent histories approach applied to systems 
consisting of large numbers of particles, with histories consisting of projections onto local 
densities, with the aim of showing, first, that these histories decohere; and second, their 
probabilities are peaked about hydrodynamic equations. In this paper we make a number 
of steps in this direction. 

Generally, we are concerned with systems which are described at the microscopic level 
by a Hamiltonian of the form 



The local densities of interest are the number density n(x, t), the momentum density g(x, t) 
and the energy density h{x.,t), defined by. 



When these objects are raised to the status of operators in the quantum theory, pairs 
of non-commuting operators are symmetrized to ensure hermiticity. (In what follows it 
will be assumed, where necessary, that non-commuting operators are ordered in this way, 
even though it will not always be written explicitly). The standard derivation of the 
hydrodynamic equations for the average values of the local densities starts by considering 
the continuity equations expressing local conservation [12]. These are of the form 




(1.1) 




(1.2) 




(1.3) 



(1.4) 




(1.5) 
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where a denotes n, g or h (and the current j is a second rank tensor in the case of g). 
It is assumed that, for a wide variety of initial states, conditions of local equilibrium are 
established after a short period of time. This means that on scales small compared to the 
overall size of the fluid, but large compared to the microscopic scale, equilibrium conditions 
are reached in each local region, characterized by a local temperature, pressure etc. which 
vary slowly from region to region. This state of affairs is described by a local equilibrium 
density operator, which has the form 

p = Z-^ exp y d^a;/?(x) [/i(x) - /i(x)n(x) - u(x) • g(x)]^ (1.6) 

where and u are Lagrange multipliers and are slowly varying functions of space and 
time. (/? is the inverse temperature, u is the average velocity field, and /i is related to the 
chemical potential which in turn is related to the average number density). The hydro- 
dynamic equations follow when the continuity equations are averaged in this state. The 
averaging produces averages of the currents, (j), but constitutive relations then typically 
follow, which relate the averages of the currents to averages of the local densities (and 
their derivatives). Hence one obtains a closed set of evolution equations for (n), (g) and 
(/i) (smeared over small spatial volumes). 

What is therefore known already is that the average values of the local densities, in 
the local equilibrium state, obey hydrodynamic equations. 

Demonstrating emergent classicality, however, requires considerably more than this. 
The above result is not unlike the Ehrenfest theorem of elementary quantum mechanics. 
This shows that the expectation values of position and momentum obey classical equations 
of motion, but it alone does not constitute a derivation of the emergence of classical 
mechanics from quantum mechanics [13]. 

Here, we adopt the discussion of emergent classicality used in the decoherent histories 
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approach. In that approach we consider the decoherence functional, 

Dia,g!) = Tt [PaM ■ ■ ■ PaAh) pPa'S^l) ' ' ' Pa'Jtn)) (1-7) 

The histories are characterized by projectors Pa{t) (in the Heisenberg picture) at each 
moment of time, in this case onto the local densities (the detailed construction is described 
below). The decoherent histories approach is described in detail in many places, so we 
will not go into detail here [1,2,3,4,5,14]. Briefly, the histories are said to be decoherent 
when the decoherence functional (1.7) is approximately zero in its off'-diagonal terms. 
Probabilities equal to the diagonal elements may then be assigned to the histories in the 
set. 

A decoherent histories derivation of the hydrodynamic equations goes beyond the stan- 
dard one in two respects. First of all, it shows that histories of local densities are decoher- 
ent, so probabilities can be assigned to them. Secondly, rather than considering just the 
evolution of average values of the local densities, we consider the probability for histories of 
values, which is a far more detailed object. We will, however, find an important simplifying 
feature in the special case of a local equilibrium initial state: this is that the probabilities 
for histories are strongly peaked about the evolution of the average values. In this case, 
the jump from the evolution of average values to probabilities for histories is then in fact 
not that large. 

It is perhaps also of interest to note that it is difficult to discuss the decoherence and 
emergent classicality of the local densities using other approaches to decoherence. The 
density matrix approach, for example, relies quite heavily on the existence of a reduced 
density matrix for the variables followed, e.g., the position of Brownian particle, after the 
variables ignored are traced out. This requires that the total Hilbert space is a tensor 
product of the Hilbert space of the variables followed with the Hilbert space of the en- 
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vironment, which is not the case here. The decoherent histories approach, by contrast, 
adapts more comfortably to these situations. It deals with projections onto subspaces of 
the total Hilbert space, and does not require a tensor product structure. 

Turning to this paper, we will address the issue of the emergence of hydrodynamic 
equations for the special case of the diffusion of local number density of foreign particles in 
a fluid. We then go on to argue that the approach and some of the results also apply to the 
more general case; thus, a (partially heuristic) derivation of the hydrodynamic equations 
from decoherent histories will be given. 

The system we are initially interested in is a fluid containing a very low concentration 
of foreign particles. The number of foreign particles is conserved, and it is well-known 
that the local number density obeys the diffusion equation. For us, the simplifying feature 
is that the microscopic dynamics of the foreign particles is reasonably simple, and the 
approach to equilibrium and the consitutive relations are easily derived. In particular, 
each foreign particle may be regarded as a free particle in a bath (the background fluid), 
so is essentially a particle undergoing quantum Brownian motion, the description of which 
is well-known. 

Hence, in this particular model there is in fact, an environment, namely the background 
fluid, which assists in many ways, in particular it makes a contribution to decoherence. 
What is new, however, compared to previously considered system-environment models, is 
to consider histories of local densities of many Brownian particles, rather than histories of 
positions of a single point particle. Indeed, this model is a kind of half-way house between 
the familiar system-environment models and true hydrodynamic models. Furthermore, as 
indicated above, the model is sufficiently similar to the proper hydrodynamic case to make 
its results of relevance. 
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The quantum Brownian motion model and its diffusive limit is described in Section 
II. Its reinterpretation as a deterministic equation for the number density is described in 
Section III. The decoherent histories analysis is described in Section IV. It is shown that 
histories are peaked about diffusive evolution, and that they are approximately decoherent. 
The extent to which these ideas go over to the more general case is discussed in Section 
V. We summarize and conclude in Section VI. An alternative proof of the key results of 
Sections IV and V is given in an appendix. 

A number of other authors have investigated closely related issues. The present work 
is strongly rooted in the decoherent histories approach to quantum theory of Gell-Mann 
and Hartle, who have frequently stressed that the most general situations in which there is 
decoherence and peaking of the probabilities about deterministic equations are generally 
not of the system-environment type [2,3,13]. An earlier paper investigated the decoherence 
of histories of number density in a simple model involving spin waves [15]. Brun and Hartle 
have investigated coarse grainings in an oscillator chain model which allow a wide variety 
of different divisions into system and environment, and also lead to a wave equation as 
the effective classical equation of motion [16]. Calzetta and Hu have investigated the 
decoherence of histories specified by fixed values of certain correlation functions in a field 
theory [17]. Other coarse graining schemes not of the system environment type (and 
not involving decoherent histories) may be found in Refs. [18,19]. Some related issues 
concerning the derivation of hydrodynamic equations may be found in Ref. [20] . Omnes [21] 
has observed that the work of Feffermann [22] may contain some useful guides concerning 
the selection of a preferred set of variables is a general system. 

Finally, it should be mentioned that not all authors agree that the decoherent his- 
tories approach is capable of explaining the emergence of classical physics from an un- 
derlying quantum theory. The various sides of this rather subtle debate are described in 
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Refs. [23,24,25]. The present paper is essentially concerned with exploring the mathemat- 
ical properties of the decoherent histories approach, as it currently exists, and has little 
to add to this debate. The extent to which the decoherent histories approach remains a 
valuable and useful one is perhaps best summarized in Ref . [26] . 



2. DERIVATION OF THE DIFFUSION EQUATION 

The diffusion equation arises for a system consisting of a collection of N foreign particles 
in a fluid of a much larger number of background particles. We assume that the density of 
foreign particles is so low that their interactions with each other may be neglected. Each 
particle interacts only with the fluid which we assume is in equilibrium at temperature T. 
Each individual foreign particle, therefore, may be modeled as a free Brownian particle for 
which the fluid plays the role of a thermal environment. The description of this system is 
well known. Classically, the phase space density W{pj q) of the particle obeys a Fokker- 
Planck equation of the Kramers type, which takes the form [27] , 

In the quantum theory, the Wigner function obeys the same equation (in the high 
temperature limit, which we work in here). The density operator is recovered from the 
Wigner function by the transformation 

p{x, y) = jdp et^'(^-y) W{p, ^) (2.2) 

and it obeys the master equation 

dp ih [ d'^ p d'^ p 



dt 2M \dx'^ 

[dp dp\ 2M7fcT 2 
- 7(- - y) - ^ J - -^(x - y) P (2.3) 
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This equation arises directly in the Caldeira-Leggett model of quantum Brownian motion, 
in which a particle is linearly coupled to a thermal bath of harmonic oscillators [7] . How- 
ever, there are many indications [28,29] that the form of this equation is in fact a lot more 
general, and will still hold, in some approximation, for more complicated and realistic 
couplings. 

Eq.(2.1) may be solved by introducing the propagator K{p, q, t\pQ, qo, 0), defined by 

W{p, q,t) = J dpodqo K{p, q, t\po, qo, 0) Wo{p, q) (2.4) 
The propagator also obeys Eq.(2.1) with the initial conditions, 

K{p, q, 0\po, qo, 0) ^ S{p - po) 5{q - qo) (2.5) 

The explicit form is, 

K{p, q, t\po, qo, 0) = exp (-a(p - p^^f - (5{q - q^^f - e{p - p^^){q - g^^)) (2.6) 
(ignoring prefactors, which may be recovered by normalization where necessary). Here, 

p"* = (2.7) 

Explicit expressions for the coefficients a, (3, e may be found in Ref.[30]. 

We are interested in the behaviour of this system for times t » Unlike the 

harmonic oscillator, this does not settle down to equilibrium, but certain simplifications 
do occur. First of all, note that 



2M7 



^ (2. 



Second, from the explicit expressions in Ref . [30] , it can be shown that 

e^-^ (29) 

The solution Eq.(2.1) therefore approaches the form, 

W{p, q,t) = J dpodqo exp (-ap^ _ _ ^d)2 _ ^^(^ _ ^d)^ ^^(^^ ^) (2.IO) 

The distribution of momenta is 

g{p) = jdq W{p, q, t) ^ exp (-^^) (2.11) 

(where we have used ^t» 1). Hence the momentum distribution approaches the Maxwell- 
Boltzmann form. The distribution of positions, on the other hand, approaches the form 

/(g,0 = j dpW{p,q,t) 

^ J dpodqo exp - go - ^) j Wq{po, qo) (2.12) 

It follows that f{q,t) obeys the diffusion equation, 

where D = kT/2M-f. 

This equation can also be derived by integrating Eq.(2.1) over p, with the result, 

and Eq.(2.10) may be then be used to show that, 

/kT d 
dppW{p,q,t)^- ——f{q) (2.15) 

hence the diffusion equation follows. 
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3. DIFFUSION OF NUMBER DENSITY 



Turn now from the case of a single particle to the system of interest - a collection of 
N particles, each one of which is described by the one-particle dynamics outlined above. 
The number density at the point x is 

N 

^(^) = ^(lk-x) (3-1) 
k=l 

Usually, we are interested in the number density smeared over a small volume V, which 
we write, 

TV 

nv{x) = Yl ^viQk - x) (3.2) 

k=l 

where Sy denotes a delta function smeared over V, i.e., a window function equal to 1 in an 
interval of size V and outside, and x will then be a discrete label. (We will not explicitly 
use the smeared form unless necessary). Microscopically, we assume that the N particle 
system is described by the phase space distribution function 

N 

/^(P1,P2,---P7V,91,92,---97V) = n ^(Pk^lk^t) (3.3) 

k=l 

where W{p, q, t) are the one-particle Wigner functions of the previous section. The prod- 
uct form is reasonable since the particles do not interact, and even if they are initially 
correlated, these correlation will eventually be lost to the environment of each particle. 

The number density n{x) is a function of the stochastic variables q]^. However, we will 
show that the evolution of n{x) is in fact strongly peaked about deterministic evolution 
according to the diffusion equation. This is, of course, almost obvious, since all we are 
doing is in essence reinterpreting the one-particle probability as a distribution function on 
a real ensemble of N particles. But for what follows, it is worth going through the steps. 
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First we show that the average value obeys the diffusion equation, then we show that the 
distribution is strongly peaked about the average. 

Consider first the average of n{x) in the N particle distribution (3.3). It is, 

{n{x)) = J2 {^{'lk-x))N 
k 

= N{5{q-x))i 

= N j dpW{p, X, t) = Nf{x) (3.4) 

where {■ ■ ■)n ^^'^ {' ' ')l denote averages in the A^-particle and one-particle distributions, 
respectively. Hence {n{x)) obeys the diffusion equation, since f{x) does, as shown in the 
previous section. 

It is, however, enlightening to briefly demonstrate this a different way, more closely 
related to standard derivations of hydrodynamics equations. The number density obeys 
the local conservation law, 

dn 1 do , , 

where 

9{x) = J2 Pk ^{ik- ^) (3-6) 
k 

Averaging through (3.5), we find we need to compute 

k 

= N J dpp W{p,q,t) (3.7) 

We have already computed this in Eq.(2.15). Hence we derive the constitutive relation, 

kT r) 

{g{x)) = -— -{n{x)) (3.8) 
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and the diffusion equation for {n{x)) follows. Note also that the state averaged in here is 
close to a local equilibrium state - the momenta are in a Maxwell-Boltzmann distribution 
(2.11), the temperature T is constant throughout the system, the average velocity flow is 
zero, and only the number density varies throughout the system (although the correlation 
between p and q in (2.10) is important for (3.8) to be obtained). 

Now consider the variance of n{x) in the distribution (3.3). Here we have to consider 
the smeared from (3.2). We have 

(ny) = J2 {^viQk - x))n = N{6v{q - x))i (3.9) 
k 

and 

jk 

= N{S^{q - x))i + (iV2 - N){5v{q - x))l (3.10) 

Therefore, 

(ny)2 -AT (5y)2 ^''■''^ 
For large N, which we assume, the distribution of ny(a;) is very strongly peaked about its 
mean value. 

The above result, which is a form of the central limit theorem, is important in the 
derivation of more general hydrodynamic equations (and indeed is used again in the next 
section in a slightly different form) . It essentially relies on two features of the macroscopic 
system. Firstly, that the microscopic distrbution is of the form (3.3), i.e., a product 
of identical one-particle distribution functions. Secondly, that the macroscopic variable 
in question, such as (3.1), is a sum of identical terms each of which involves only the 
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one-particle variables. It is, however, very likely (and probably provable) that the result 
will hold more generally. In particular, we expect it to continue to hold in the more 
realistic situation where the microscopic distribution function involves a certain amount 
of correlation between particles (so is not exactly of the product form) and when the 
macroscopic variables involve interaction terms between particles. Generally, we expect 
strong peaking about the mean values of a macroscopic variable when fixed values of that 
variable correspond to a very large number of similar microstates. 



4. HISTORIES OF NUMBER DENSITIES 



4(A). Probabilities and Decoherence 



Consider now the question of histories of number densities. After tracing out the envi- 
ronment for each of the N particles, the decoherence functional for histories characterized 
by projections onto number density at two moments of time is 



(4.1) 



The trace is over the Hilbert space of the N particles. Kq is a product of N reduced density 
matrix propagators (corresponding to the master equation (2.3)), one for each Brownian 
particle. Pn denotes a projector onto number density, and may be represented as 



\x) = Hx) - n{x)] 



(4.2) 



Here, So- is a delta function smeared over a range a. The number density operator n{x) may 
be smeared over a range V when necessary. (More detailed representations of projections 
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onto number density are possible, in terms of projections onto position of the individual 
particls [15], but this is not necessary here). Eq.(4.1) is the decoherence functional for the 
simplest non-trivial case of two-times histories, but the arguments present below will also 
go through for histories of projectors at many moments of time. 

For the initial state p, let us again suppose that it is a product of N identical terms, 
as in the previous section. We will also assume that it has already evolved for a time 
t >> 7~^, so we are in the diffusive regime. We will consider later the consequences of 
other types of initial state. 

If the initial density matrix is a product of N identical terms, we expect a central limit 
theorem result to hold, and the distribution of n{x) in the initial state will be strongly 
peaked about its mean. In more detail, suppose we have a pure initial state on the N 
particle Hilbert space which is a product of N identical terms, 

|*) = |V')«)|V')«'---«'|V') (4.3) 

Then for large N, \^) is approximately an eigenstate of the number density operator (2.1), 

n{x)\^) ^ {nix))\^) (4.4) 

where {n{x)) = N\'i/j{x)\'^. This is essentially the Finkelstein-Graham-Hartle theorem [31], 
since the number density is the same as the relative frequency operator (up to an overall 
factor of N). For projections onto number density a similar result follows, that is, 

6a [fi{x) - n{x)] I*) fa 6a [{n{x)) - n{x)] [*) (4.5) 

It is straightforward to generalise to the case of mixed states. That is, for a density matrix 

p = Pl ® pi ® ■ ■ • ® pi (4.6) 
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it may be shown that 

n{x)p ~ {n{x))p (4.7) 
and thus, the analogue of (4.5) holds. 

These results mean two things. First, that 

Pn,pPn'^-0, for ni ^ n[ (4.8) 

Hence, there is approximate decoherence for initial states of this type. Secondly, they allow 
us to compute the probabilities. We have 

PnipPfii ~ Sa [(n(a;)) - ni{x)] p (4.9) 
Inserting into the diagonal elements of the decoherence functional (4.1), we find 

p{n2, ni) = Tr {Pn.pt) 5^ [{h{x)) - ni{x)] (4.10) 

Now pt is the evolution of the initial state under the non-unitary propagator Kq. The key 
point is that it remains in the product form, so, 

TV {Pn,pt) ^ Sa mx))t - Mx)] (4.11) 

Hence, the probability (4.10) is strongly peaked about ni{x) = {n{x)) and n2{x) = (n(a;))f, 
that is, about evolution according to the diffusion equation. 

We have therefore proved the desired result: for initial states of the form (4.6), histo- 
ries of number density are decoherent, and their probabilities are peaked about diffusive 
evolution. 

4(B). More about Decoherence 

Let us consider now in more detail the question of decoherence. Since each individual 
foreign particle is coupled to its surrounding fluid environment, the possibility exists that 
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this environment causes decoherence of each particle's position. Decoherence of number 
density then follows trivially, since projections onto number density may be constructed 
from coarse-grainings of projections onto position [15]. However, decoherence of a particle 
coupled to an environment requires certain choices of parameters for the model, such 
as large mass and large temperature. Hence, there will be certain parameter ranges (in 
particular, small masses for the foreign particles) for which the decoherence at the one 
particle level may not be very good. In this case, it is the large N effect which produces 
the dominant contribution to decoherence of number density. (Generally, replicating the 
system N times enhances the degree of decoherence exponentially in N [15]). 

Because decoherence is not necessarily produced by the environment of each parti- 
cle, it is of interest to see how decoherence comes about for initial states consisting of 
superpositions of number density. So let, 

v2 

+ |x)«)|x)®---«'lx)) (4.12) 

It is now the case that Pfj^pP^/^ 7^ if we take ni{x) = N\ij;{x)\^ and n'^{x) = N\x{x)\^. 

In this case, decoherence may be demonstrated by appealing to approximate determin- 
ism: decoherence follows if the probabilities are strongly peaked about a unique relationship 
between alternatives at different times. To see this, note that for histories consisting of 
alternatives at two moments of time, the decoherence functional satisfies [5] 

\D{ni,n2\ni,n2)\ < p{ni,n2)p{ni,n2) (4.13) 

So if p(ni, 77,2) is strongly peaked about an equation linking ni and 712, only one of p(ni, n2) 
and p{ni,n2) can be non-zero for fii ^ n'^, hence the decoherence functional is approxi- 
mately diagonal. 
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Before fooking in more detail at how this argument apphes to the specific case of 
diffusive evolution, some qualifications must be added. For diffusive evolution there is not, 
in fact, a unique initial density corresponding to a fixed final one. There will generally be 
a spread of initial densities that evolve into a given final density, and the non-uniqueness 
will get worse with increasing time intervals. A fixed final density will fix the initial density 
precisely only for reasonably short times, and for sufficiently large choices of the width cr 
of the projections onto number density (which ensures that different values of n{x) are 
sufficiently "far apart"). Fortunately, these restrictions do not adversely affect the main 
objective here - to demonstrate the emergence of the diffusion equation - since it is clearly 
sufficient to consider projectors onto number density at two closely separated times. 

Returning to the demonstration of decoherence, the probabilities pini,n2) are non-zero 
only if ni is equal to either N\ip{x)\'^ or N\x{x)\'^. The first choice projects the initial state 
onto lip) <Si \ip) <Si ■ ■ ■ and the second onto |x) <^ Ix) <^ ' • The point is now that each of these 
two states, on evolution for a short time t, are strongly peaked about two distinct values 
of number density, as determined via Eq.(2.12). (It is not necessarily diffusive evolution, 
except for times t » It is only necessary that the final value of number density 

be uniquely determined by the initial state). It follows that, modulo the above remarks, 
p{ni,n2) is strongly peaked about a unique relation between ni and 77-2, hence there is 
approximate decoherence. 

4(C). Other Variables 

It is possible to consider histories of other types of variables in this diffusive model. 
Momentum and energy are not so interesting because they are not locally conserved - the 
momentum and energy of the foreign particles is exchanged with the environment. But 
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another finer-grained variable that is hkely to be of more general use is the phase space 
density. Classically, this is given by 

fix, k) = J2 KVj - k)5{qj - x) (4.14) 
3 

This quantity may be smeared over ranges of p and q, so that it becomes a function /r(ic, k) 
which counts the number of particles in the phase space cell of size F labeled by discrete 
lables x,k. Classically, 

{f{x,k))N = NW{k,x) (4.15) 

where W{k,x) is the one-particle phase space distribution function. Like the number 
density, the distribution of f{x,k) is strongly peaked about its mean, (4.15), since (4.14) 
is a sum of identical terms. Hence, we can regard the phase space density (4.14) as a 
hydrodynamic variable whose evolution equation is the Fokker-Planck equation (2.1). 

Similar results will hold in the quantum theory. Operator ordering is necessary to 
make (4.14) hermitian. We can the consider projections onto ranges of its values. 



5a 



fr{x,k)-f{x,k) (4.16) 



The subsequent treatment will be similar to the number density case. One would expect 
the quantum theory to exhibit an approximate determinism in the evolution of (4.14) 
(similar to Omnes theorems on the evolution of phase space projectors [4]) and this will 
guarantee approximate decoherence of histories of phase space density. The probabilities 
for these histories will then be peaked about evolution according to the Fokker-Planck 
equation. (This will be investigated in more detail elsewhere). 



5. General Hydrodynamic Histories 



Interestingly, many of the features of the diffusive model go over to the general case. 
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Consider a closed system consisting of N weakly interacting particles, now without a back- 
ground fluid. The number, momentum, and energy density are given by 

^i^) ^^^{Qj (5-1) 
j 

9{x) = ^S{qj - x)pj (5.2) 
j 

2 

Kx) = Y,S{qj-x)^ (5.3) 
j 

They are written for the case of a one-dimensional model for notational simplicity, al- 
though the three-dimensional case is readily obtained. The Hamiltonian is for a system of 
approximately non-interacting particles. This is sufficient to cover, for example, the case 
of a dilute gas. The neglect of interactions might seem restrictive, but, as in the diffusive 
model, we find that the interesting things we are able to say about the probabilities histo- 
ries arise from large N statistics, rather than from detailed knowledge of the microscopic 
interactions. 

As in the previous case, the densities need to be smeared over a small volume. Also, 
in the quantum case, non-commuting operators should be ordered in such a way that they 
are hermitian. 



5(A). Decoherence and Probabilities 
for a Local Equilibrium Initial State 



The usual derivation of the hydrodynamic equations [12] involves assuming a local 
equilibrium state, which for this one-dimensional case has the form 

p = exp ^— J dxP{x) h{x) — p,{x)n{x) — u{x)g{x)^ (5-4) 
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Because of the delta-functions in (5.1)-(5.3) this becomes, 



p = Z ^ exp 




2m 



— - fiiqj) - u{qj)pj 



(5.5) 



The usual manipulations on this object may be performed by assuming that the Lagrange 
multipliers /3, /i, u are slowly varying functions. In particular, one may show that the 
corresponding one-particle Wigner function is 

Wife,) = /(,) exp (-<|I^) (5.6) 

where kT = 1/(3. The quantities f{q)iU and T are assumed to be slowly varying functions 
of q and t. (The quantity f{q) is proportional to e^^^, where fi = — ^mu^). It is 
straightforwardly shown that the averages in this state, {n{x)), {g{x)), {h{x)) obey a 
simple set of hydrodynamic equations. 

In the quantum theory, projectors onto local densities are again constructed using delta- 
functions. The projection onto number density is as before Eq.(4.2) and the projectors 
onto momentum and energy density are 

P-g = 5a, [9{x) - -g{x)] (5.7) 
h{x) - h{x)\ (5.8) 



P-h = 5<^H 

Again S^r denotes a delta function smeared over a range cr. Strictly, there is one such 
delta-function for each value of x (which is discrete if 5{qj — x) is smeared over a small 
range). Although the operators h{x), g{x), h{x) do not commute, the projections onto 
ranges of their values Pn, Pg, Pj^ will approximately commute with each other if both the 
width (7 of the projections and the smearing volume of 5{qj — x) is sufficiently large. (The 
extent to which this approximation is good will also depend to some extent on the initial 
state they are traced with in the decoherence functional). 
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In the approximation that these projectors commute, the decoherence functional for a 
two time hydrodynamic history are 

Dial, a2\a', a^) = Tv (^Pn,P-g,Pj,^e-^^* PnAPn, P Pn.Pg.Ph/"') (5.9) 

We take the initial state p to be the local equilibrium state (5.4). From what has gone 
before, it should be clear that results similar to (4.8), (4.9) in the diffusive model hold. 
That is, if the widths of the projections an, erg, (Tj^ are much greater than the variances of 
the corresponding operators in the local equilibrium state, we have 

Pn.Pg^P-h, pPn'P-g'P-h[ ^ (5-10) 
for ni ^ n'^, gi ^ g'-^, hi ^ h'^. Hence there is decoherence for this initial state. Similarly, 

PniP-giP-h, pPniP-giP-h, ~ ^'^n mx))-ni{x)] 

X ^a, mx)) - gi{x)\ 6a, [{h{x)) - hi{x)] p (5.10) 

These two results follow from the fact that the local equilibrium state is very strongly 
peaked in the locally conserved quantities about their average values for large N. The 
proof of this property of the equilibrium state is similar to the proof, in Section 3, that 
the distribution of number density in the diffusive model is strongly peaked about its 
mean. That proof relied on the property that the averages of ny and riy are sums of 
identical terms involving averages in the one-particle distribution functions (see Eqs.(3.9)- 
(3.11)). In the local equilibrium state, the averages of functions of the local densities are 
not exactly sums of identical terms, since the parameters of the distribution, such as are 
not constant. That is, the one-particle distribution functions such as (5.6) are not quite 
the same from one particle to another. The required peaking about the mean will still 
hold, however, because the parameters of the local equilibrium state are slowly varying 
functions. 
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After being hit by the initial projection, the density matrix is evolved to time t, where 
it is still a local equilibrium state, and hitting it again with a projection yields the average 
values at time t. Hence, 



n 



[{n{x))t - n2{x)] 5ag[{g{x))t- g2ix)] Saf, {h{x))t - h2{x) 



xSa„[{h{x)) -ni{x)] Sag[{g{x)) - gi{x)] Sa^ {h{x)) - hi{x) 



(5.12) 



The probability distribution is therefore peaked about the hydrodynamic equations. (An 
alternative derivation of this result, using Gaussian projections, is given in Appendix A). 

This is the main result of this section: for a local equilibrium initial condition, histories 
of local densities are approximately decoherent and their probabilities are peaked about 
evolution according to the hydrodynamic equations. 

The discussion of finer-grained histories involving projections onto phase space density, 
Eqs.(4.14)-(4.16) also applies to the more general case. Here, the phase space density is 
not unlike a locally conserved quantity, since it can change only as a result of a net flow in 
phase space over the boundary of the local phase space cell. (This is not the case in the 
diff'usive model). The main difference is that in the general case we would expect the one- 
particle phase space density to obey a Boltzmann equation, rather than a Fokker-Planck 
equation (2.1). 



Turn now to the question of decoherence for more general initial states. As stated in 
the Introduction, we expect the approximate conservation of the local densities to ensure 
their approximate decoherence. This is because exactly conserved quantities are exactly 



5(B). Decoherence More Generally 
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decoherent [11], which may be seen as follows. Consider the two-time decoherence func- 
tional: 



Suppose the projections are onto some exactly conserved quantity. Then Pq!2 commutes 
with iy, so it may be moved through the unitary evolution operator to act directly on 
or Pq,/^. Since Pa2Pai = unless 0:2 — Q^i, we see that the decoherence functional is zero 
unless ai = a2 = a'^, hence there is decoherence, for every initial state. 

More generally, there will also be (approximate) decoherence if the quantum theory 
exhibits some kind of approximate determinism. Suppose that evolution of the projectors is 
such that e^^^ P(X2G~^^^ ~ ^a2{t)- That is, it is approximately equal to another projector 
of the same operator (but onto a different range, labeled by 0:2 (t)). Then again the 
decoherence functional is approximately zero unless ai = a2{t) = a'^. An example of this 
type of decoherence may be found in the histories of phase space projectors, considered by 
Omnes. He showed that a projector onto a sufficiently large cell Fq of phase space evolves 
under unitary evolution into another phase space projector onto the phase space cell Ft, 
the classical evolution of Fq. 

In a similar way, we can argue for the decoherence of local densities for a wide class of 
initial states by considering the time evolution of the final projector. Take it to be of the 
form 



and similarly for the initial projectors, Pai- The question of decoherence is the question 
of the general form of the object. 




(5.13) 




(5.14) 



-iHt 



(5.15) 
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In particular, if this object is strongly localized in n, g and then there will be approximate 
decoherence, since 

Pa,^Pa2{t)Pa^^0, for ai ^ a[ (5.16) 

This is of course an operator statement, so involves some kind of restriction on the class 
of initial states p it is traced with in the decoherence functional. But we would expect it 
to be true for a wide class of states, although clearly not all. 

To see why (5.16) should hold, regard the final projector as representing a final state, 
Pf — Pa2/T^{Pa2) (the trace is finite). So Eq.(5.15) is the backwards evolution of this 
state, 

pj{-t) = e*%je-^^* (5.17) 

Now, for unitary evolution, evolution backwards in time has the same generic physical 
features as evolution forwards in time. In particular, it is strongly believed that a wide 
class of initial conditions for a large many particle system will quickly evolve to a local 
equilibrium state [32]. Therefore, we expect that for sufficiently large t (5.17) will be close 
to a local equilibrium state, hence will be localized in n, g and h, and so (5.16) will hold. 

The issue is further helped along by the observation that Pq,2 hence P/ is initially 
localized in n, g and h, and, since these are locally conserved quantities, both they and 
their fluctuations will vary slowly. In fact, p/ is closely related to a local equilibrium 
state. It is in a sense the microcanonical version of it. In pj, the variables n, g and h are 
constrained to lie in a small range centred about n, g and h. In the local equilibrium state, 
on the other hand, only the averages of n, g and h are constrained to take set values, but 
large N statistics guarantees that the fluctuations about the mean are extremely small, 
hence the distributions of n, g and h are here also concentrated in a small range centred 
about the mean. 
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Since p/ is either essentially equivalent to a local equilibrium state, or very close to one, 
and since a local equilibrium state remains so under unitary evolution, it is very plausible 
that (5.17) is localized in the hydrodynamic variables (and in particular it is probably a 
local equilibrium state), so (5.16) will hold. We therefore expect decoherence for a wide 
variety of initial states. 

It is perhaps pertinent to elaborate a little on what it means for a state to approach 
the local equilibrium form (5.4). We expect, for instance, a physical system starting out 
in a pure initial state to approach local equilibrium, in some sense. Yet (5.4) is a very 
mixed state so cannot be reached from a pure state by unitary evolution. What is going 
on here is that one is typically interested only in the averages of the local densities, and 
for these purposes, the quantum state of the system (which may be pure) can be replaced 
with another state possessing the same average values. Usually, the replacement state is 
chosen by the maximum entropy method (subject to fixed average values), and it is this 
which leads to the local equilibrium state (5.4) [33]. 

Differently viewed, since the local densities are sums of one-particle operators, their 
averages involve only the one-particle distributions functions, such as (5.6). (If interactions 
are significant two-particle operators and distribution functions could be involved) . These 
distribution functions may exhibit the local equilibrium form (5.6), even though the actual 
state of the whole system is not of this type. The local equilibrium density operator of the 
whole system, (5.4), is the effective state of the whole system given that the one-particle 
distribution functions are of the form (5.6). 

In this connection, it is also worth mentioning a different argument for decoherence. 
We are interested in histories of projections onto number, momentum and energy density 
on the system of N particles. However, as indicated in Ref.[15], these projectors may 
be expressed in terms of a (generally rather complicated) sum of projections involving 
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each individual particle. Hence, we could compute the decoherence functional for the local 
densities by first considering the decoherence functional for the much finer-grained histories 
involving the position, momentum density and energy density of each particle. Now the 
point is that the decoherence functional for each individual particle would then have the 
familiar system-environment form, in which the properties of one or two particles are 
followed and all the rest are traced out. On general grounds, we would therefore expect a 
degree of decoherence of the properties of each individual particle. Since the local densities 
of the N particle system are constructed from coarse grainings of the individual particle 
histories, the degree of decoherence can only improve (as seen explicitly in Ref.[15]). From 
this point of view, the hydrodynamic models are in fact very similar to the difi^usive model 
of this paper, except that here each particle doubles up as both system and environment. 

The main difference, however, between the hydrodynamic models and the usual system- 
environment models, is that there are no distinguished particles, since they will typically all 
have the same mass. So although there may be a degree of decoherence in the one-particle 
reduced density matrices, unlike the quantum Brownian motion models, there will be little 
determinism in the evolution of the single particle positions and momenta, because the 
fluctuations will be very large as a result of interacting with an "environment" consisting 
of large numbers of particle of the same mass as the "system". This, in fact, is one of 
the reasons it is interesting of consider instead histories of local densities of the N particle 
system - averaging over very large numbers of particles smooths out the fluctuations, so 

that the local densities evolve almost deterministically. 
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6. SUMMARY AND CONCLUSIONS 



For a local equilibrium initial state, the derivation of hydrodynamic equations from 
decoherent histories is reasonably simple. Large statistics guarantees that the prob- 
abilities for histories of local densities are strongly peaked about a single history - the 
evolution of the mean in the local equilibrium initial state. We therefore readily connect 
with previous work, which shows that the mean values of local densities evolve according 
to the hydrodynamic equations. Furthermore, the fact that there is essentially only one 
history with non-zero probabilty means that there is approximate decoherence. 

Going beyond the assumption of a local equilibrium state, is however, much more 
difficult. We were able to see how this sort of state arises asymptotically in the simple 
diffusive model, but it is difficult to demonstrate in general. We gave some arguments for 
decoherence for more general initial states, but these still remain rather heuristic, and it 
would be useful to develop some explicit models in which decoherence via approximate 
conservation could be seen more clearly. This will be addressed in future publications. 
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THEOREM FOR HISTORIES 



In this appendix, we give an alternative derivation of the result of Section V, that the 
probabilities for histories of local densities are peaked about the evolution of the mean 
values (and hence about the hydrodynamic equations). There, the result relied on large 
N statistics. Here, we show that it may be regarded as an example of a simple theorem 
for histories similar to the Ehrenfest theorem of elementary quantum mechanics. The 
result is as follows: for any Hamiltonian and any initial state, histories characterized 
by projections onto any operators are decoherent if the widths of the projections are 
taken to be sufficiently large, and the centers of the projections are also suitably chosen. 
Furthermore, the probabilities for these histories are then peaked about the time evolution 
of the expectation value of the operators in the given initial state. 

At some level, this result must be almost trivial, since it is so very general ( "any" initial 
state, "any" Hamiltonian, "any" operator). However, as we shall see, the hydrodynamic 
case considered in this paper is a situation to which it naturally applies, where it is little 
more than the statement that the quantum evolution is strongly concentrated around one 
particular history. 

To give a flavour of the basic idea, consider the simplest possible case of a history 
characterized by a single projection at one moment of time. Such histories are trivially 
decoherent, and the probabity is given by the standard expression. 



For simplicity, we will let the projection Pa be a Gaussian quasi-projector onto some 
operator A with width a: 




(Al) 
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Gaussian projectors are exhaustive but only approximately exclusive. 

Inserting this in (Al) and letting P = a — (A), and B = A — (A), we have 



1 , ( iP-B 



We are interested in the behaviour of this expression for large a. Expanding for large a 
we obtain, 

(note that P is taken to be of order cr because of the Gaussian factor). Now, noting that 
(B) = 0, and = {AAf, and reorganizing the result as an exponential, we find that 
for large cr. 

Pi-) - ' . . exp (^5) 



(27r(a2+(AA)2))i V 2(a2 + (AA)2), 
Therefore, for » (A^)2, the probability is peaked about the expectation value of 
A. This is an intuitively expected result: a very broad sampling yields the average value. 
(The same qualitative result will probably also hold for exact projectors, but we do not 
consider this here). 

The above result readily extends to non-trivial histories. First note that for Gaussian 
projectors, the time-dependent projectors may be written. 



(27ra2) 

where 



Then inserting in the expression for the probabilities (setting aside for the moment the 
question of decoherence) and expanding for large cr, the leading order behaviour is, 

p{ai,a2,---an) = fl exp {A8) 
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For simplicity, we have ignored normalization factors. Also, we have kept only the leading 
order terms for cr^ >> (AA^j,)^ (for times t^). The next correction may be computed, at 
some length, but the qualitative conclusion is the same. 

Eq.(A8) indicates that the probabilites for histories are peaked about the evolution of 
the expectation value, {Af). It is the natural generalization of (A5): the broadest possible 
sampling of a history, yields the evolution of the expectation value. 

The result for the single time case extends easily to the multi-time case even though 
the Af^^s do not commute, because their non-commutativity only enters at higher orders 
for large u, and does not affect the configuration {Atf,) about which the probabilities 
are peaked. For the same reason, one can also project onto a number of different non- 
commuting operators at each moment of time, as one needs to in the hydrodynamic case 
(see (5.7)-(5.9)), and still obtain a result of the form (A8). 

The argument for decoherence is similar to those in Sections IV and V - it arises 
because there is essentially only one history with non-zero probability. More precisely, 
define a set containing just two histories as follows. Take the first history to be that 
characterized by projections in which the a/j's are chosen to be very close to (^i^,), and 
the width is much greater than [AA)'^. Then the probability for his history, which 
we denote simply p, will be very close to 1. Then take the remaining history to be the 
complement of the first one. That is, the combination of all histories defined by projections 
in which at least one of the ctj^'s is much greater than the width a away from (^t^). Then 
the probabilities for the second history, denote it p, will be close to zero, since at least one 
of the projections will always pick up an exponentially small tail. There will therefore be 
approximate consistency, since the probabilites approximately add up to 1. (Equivalently, 
the decoherence functional is very small, since |-Dp < pp.) 
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As stated at the beginning, the above result is very general to the point of being 
almost vacuous. But it is useful in the hydrodynamic case because the local densities do 
have extremely small variances (AA)^ (compared to (A)-^) in the local equilibrium state, 
so the widths a do not need to be taken to be unreasonably large for the conditions of the 
theorem to be satisfied. Hence we have provided an alternative proof that the probabilities 
for histories of local densities are peaked about the hydrodynamic equations. 
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